function [] = test_ts1_ts2_ts3()

amiga_model = 'a1200';
fs = get_sampling_fq();
spectrum_cut = get_spectrum_cut();

[C1, R1, C2, R2, C3 ,R3] = get_amiga_filter_params(amiga_model);

ts1 = tf(1, [R1*C1 1]);
ts2 = tf(1, [R2*R3*C3*C2 C3*(R2+R3) 1]);
ts3 = ts1*ts2;

[k, t] = firwin('a500', fs, spectrum_cut, 2048);

k=k';
t=t';

k1 = lsim(ts1, k, t);
k2 = lsim(ts2, k1, t);
k3 = lsim(ts3, k, t);

hold on;
plot(k2, '.');
plot(k3, 'r');
legend('ts1,ts2', 'ts1*ts2');

error_calc(k2', k3');

end





